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Abstract. We have measured the p t -integrated quadrupole component of two-particle azimuth correlations 
(related to quantity V2, denoted in this case by V2{2D}) via two-dimensional (2D) angular autocorrelations 
on (rj, <f>) for unidentified hadrons in Au-Au collisions at 62 and 200 GeV. The 2D autocorrelation provides 
a method to remove non-quadrupole contributions to «2 (conventionally termed "nonflow") under the 
assumption that such processes produce significant dependence on pair-wise relative rj within the detec- 
tor acceptance. We hypothesize, based on empirical observations, that non-quadrupole contributions are 
dominated by minijets or minimum-bias jets. Using the optical Glauber eccentricity model for initial-state 
geometry we find simple and accurate universal energy and centrality trends for the quadrupole compo- 
nent. Centrality trends are determined only by the initial state (impact parameter b and center-of-mass 
energy -^/ijvjv )■ There is no apparent dependence on evolving system dynamics (e.g., equation of state 
or number of secondary collisions) . Our measurements of the quadrupole and non-quadrupole components 
have implications for the contributions to v%. They suggest that the main source of the difference between 
i>2{2} and i>2\4} (or «2{2D}) is measured properties of minijets. 
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1 Introduction 



The azimuthal quadrupole in the distribution of particles 
in heavy-ion collisions, or "elliptic flow" has long been 
described in terms of hydrodynamical phenomena [lj. In- 
terpreting traditional elliptic flow measurements is diffi- 
cult because of issues with "nonflow" , particle correla- 
tions caused by phenomena such as minijets, resonances, 
and HBT. We present a two-dimensional autocorrelation 
analysis in which the quadrupole (related to elliptic flow) 
and non-quadrupole (caused by non-flow) components of 
two-particle correlations are measured simultaneously. 

In this analysis we use the structure of nonflow corre- 
lations on pair- wise relative pseudorapidity to distinguish 
geometrically the azimuth quadrupole moment, attributed 
to elliptic flow, from nonflow effects, which we believe to 
be due primarily to minijet correlations — fragments from 
low-Q 2 partons, mainly gluons [4, ; 5,,6][7llB]- Quadrupole 
amplitudes are obtained for the first time from fits to 2D 
angular autocorrelations on azimuth and pseudorapidity. 
Measurements of the quadrupole moment over a broad 
range of centralities and energies provide qualitatively new 
insights into the phenomenon of elliptic flow. Quadrupole 
amplitudes follow simple systematic trends on centrality 
and energy described by only two initial-state parameters 
for all systems down to 13 GeV. 



2 Azimuth Autocorrelations 

We construct minimum-bias angular autocorrelations by 
considering all possible pairs of particles in an event. In 
general one can consider structures in the pair density on 
6D momentum space (pti, r)i, <pi,Pt2, V2, 4>2)- In this anal- 
ysis we study p t -integrated correlations on the angular 
subspace (771, 4>i, 772, ^2), where the angle parameters for 
relativistic collisions are pseudorapidity 77 (related to polar 
angle 9) and azimuth <p. In this section we will concentrate 
on autocorrelations on only azimuth. 

For the event-wise azimuth density p(<f) the event-wise 
two-particle density is defined as P2{4>i, <fe) = p{4>i) p(4>2) ■ 
The azimuth autocorrelation is then just the projection by 
averaging of the two-particle density onto the difference 
axis 4>a = 4>i — 4*2 [2] . Consider a binned system with bin 
width S. The autocorrelation is then 

1 rit r<t>A+&/2 

r P {<t>A) = T1 #W WapWJpM) (1) 

1 r4>A+S/2 ^ ,^p± 

If we use the periodicity of 4> then in the limit where the 
bin width goes to zero we get 

r p {4> A ) = ±- p(<t>')p(<t>' ~ Mdcf>', (2) 
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which is the standard definition of the autocorrelation. 

While the autocorrelation procedure is mathematically 
well-defined for a single event, in practice the number of 
particles in an event is far too small — even in central Au- 
Au collisions — to measure two-particle densities with any 
statistical accuracy. It is therefore essential to construct 
an average of the autocorrelations over many events. The 
stationarity property of two-particle correlation structures 
[2] ensures that even in an ensemble average the informa- 
tion on the difference axis of a two-particle density — the 
autocorrelation — is not lost. The key feature of using au- 
tocorrelations is that all of the structure in an event is 
preserved, which includes both flow and nonflow effects. 
We will discuss how to seperate the two in the next sec- 
tion. 

Contrast this with the standard (event plane) tech- 
nique of measuring v 2 [3] which depends on shifting the 
single-particle densities according to the estimate of the 
reaction plane and averaging the shifted single-particle 
densities. If the reaction plane is known perfectly then 
this is an effective approach as the nonflow should not 
be correlated with the reaction plane and when averaged 
over many events will contribute as a flat background. 
However, when the reaction plane is estimated using par- 
ticles in the event then there are correlations between the 
particles used to estimate the reaction plane and the par- 
ticles used to measure the azimuthal anisotropy making 
the event-plane technique effectively another type of two- 
particle correlation |3] . This is the origin of nonflow con- 
tributions in the event-plane method. 

A proper statistical measure helps us understand how 
these correlations scale. We will refer to the pair density as 
p and define p re f to be a reference pair density constructed 
from mixed events. We can construct a per-pair density ra- 
tio Apj = p/p Te f— 1 which varies as 1/ n c h with chang- 
ing system size in the absence of other physical changes. 
We prefer a statistical measure whose variation reflects de- 
viations from the null hypothesis based on N-N superpo- 
sition, the per-particle density ratio Apj ^/p re f (Pearson's 
normalized covariance coverted to a density ratio) which 
exhibits the desired properties, since ^fp~^ oc n c h- 



3 2D Autocorrelations 

We now consider the more general problem of angular cor- 
relations on {rji,rj2, 4>\,4> 2 ). By following the procedure de- 
scribed in the previous section simultaneously on <f> and rj 
we can construct 2D angular autocorrelations on (774, 4>a)- 
Examples from 62 and 200 GeV minimum-bias Au-Au col- 
lisions are shown in Fig. Q] 

Pi-integrated 2D angular autocorrelations contain two 
types of structure: 77^-dependent peaks and ^-indepen- 
dent sinusoids cos(^>a) and cos(24>a), where the cos(20/i) 
sinusoid is related to elliptic flow, and the m = 2 Fourier 
component of the remaining structure is the source of non- 
flow in a standard two-particle analysis. We therefore de- 
compose the density ratio into a part varying with t\a 
(nonflow, identified with subscript nf ) and a part indepen- 
dent of tja composed of the m = 1 , 2 terms in a Fourier 



series 




A P _ A Pnf 



cos(m</>/i) (3) 



The correlation data for each energy and centrality 
were fitted with a six-component model function moti- 
vated by structures observed in data from 200 GeV p- 
p collisions [I]. The only addition necessary to describe 
heavy ion data is the quadrupole component. 

The model function includes: a same-side 2D Gaussian 
on (t)a, 4>a), and ry^-independent negative dipole — cos(c^a) 
an 77/i-independent quadrupole cos(20,/\) (identified with 
elliptic flow [9]), a 0^-independent ID Gaussian on T]a 
(related to participant-nucleon fragmentation), a narrow 
same-side 2D exponential on {r]Ai 4>a) (modeling quantum 
correlations and conversion electron pairs), and a constant 
normalization offset. The model function is expressed as 



_i (ha y 

F = A cf)A cos{(/) A ) + A 2<M cos(2(f> A ) + A a e n«»J 



A 2 e 



A, 



(4) 



4 Quadrupole Centrality and Energy 
Systematics 

This analysis is based on 6.7M and 1.2M Au-Au collisions 
at y/SNN = 62.4 and 200 GeV respectively, observed with 
the STAR time projection chamber [13]. The acceptance 
was defined by transverse momentum p t > 0.15 GeV/c, 
1 17 1 < 1 and 2-7T azimuth. Au-Au collision centrality was 
defined as in 12J. Minimum-bias event samples were di- 
vided into 11 centrality bins: nine ~ 10% bins from 100% 
to 10%, the last 10% divided into two - 5% bins. The 
corrected centrality of each bin as modified by tracking 
and event vertex inefficiences was determined by relating 
the mean nj^ of the bin to path length v according to 
running-integral procedures in [12j . 

The quadrupole term of our fits can be related to the 
conventional v 2 measure according to [5] 



Ap[2] 



nv\{2D} 
2~k ' 



(5) 



Fig.[2]summarizes Eq. ([4} 2D fit results for Ap[2]/ «Jp^£, 
with corresponding values of v 2 {2D} in Fig.[3]for compar- 
ison to previous analyses. The left panel shows fit results 
for 200 GeV (solid dots) and 62 GeV (solid upright trian- 
gles) data: strong increase with centrality to mid-central 
collisions followed by reduction to zero for central colli- 
sions. «2{EP} data from NA49 (inverted solid triangles) 
|14j provide a reference for energy-dependence systemat- 
ics. 

Published data for two-particle «2{2} (open circles) 
and four-particle cumulant i>2{4} (open squares) at 200 
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Fig. 1. Perspective views of two-dimensional charge-independent correlations Ap/^/p le ^ on (j)a,4>a) for Au-Au collisions at 
y / ijvjv = 200 and 62 GeV (upper and lower rows respectively) . Centrality increases left-to- right from most-peripheral to most- 
central. Corrected total cross-section fractions are noted for each centrality bin [12] . 
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Fig. 2. Quadrupole component Ap[2]/^p Te f = X[2] of the 
autocorrelation density ratio inferred from model fits to 2D 
angular autocorrelations as in Fig. [1] (solid dots and upright 
triangles) vs. centrality v (mean participant path length). 



GeV [15] are compared to i>2{2£>} (solid points) from 
this 2D autocorrelation analysis. The U2{l-D} (open tri- 
angles) are fits of cos(2^) to ID projections onto 4>a 
of the 200 GeV 2D autocorrelations, roughly consistent 
with the «2{2} analysis as expected [9lfT0] and substan- 
tially larger than the 2D fits. That nonflow offset is ex- 
pected in a conventional ID flow analysis: the difference 



rsO.09 

> 0.08 
0.07 
0.06 
0.05 
0.04 
0.03 
0.02 
0.01 





N-N STAR Preliminary 



II L 



J I L 



J I L 



1 



Fig. 3. 
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V2{2D} from 2D autocorrelations (solid points). V2 
17 GeV are inverted solid triangles in both 

Dashed 



data at ^/sjvjv 

panels. The open symbols are explained in the text, 
curves have the same form for all energies above 13 GeV (see 
text). Dash-dot curves include "nonflow." The vertical solid 
lines at v ~ 5.95 indicate v for 6 = 0. The dash-dot line at 
v — 1.25 corresponds to N-N collisions [12] . 



between open triangles and solid dots is predominantly 
the m = 2 Fourier component of the same-side minijet 
peak 10J . i>2{4} is expected to eliminate nonflow assuming 
that elliptic flow is a collective property of many particles, 
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whereas "nonflow" describes independent "clusters" of a 
few particles pj]. The open squares in Fig. [2] are closer to 
the 2D analysis, but systematic deviations outside pub- 
lished uncertainties remain. The difference «2{2} — V2{4} 
has been explained by a combination of nonflow and flow 
fluctuation effects [16J . However, we see a similar differ- 
ence in i>2{l-D} — V2{2D}, which are both two-particle 
correlations and respond identically to flow fluctuations. 
We then conclude that the main source of the differences 
in these quantities is the same-side minijet peak, not V2 
fluctuations. 

The comparison of V2{2D} and V2{1D} allows us to 
calculate the contribution of nonflow to V2{1D} in 200 
GeV Au-Au collisions to be between 20-f00%, depending 
on centrality. The problem is largest in the most central 
collisions, since not only is V2{2D} going to zero but the 
amplitude of the nonflow term is at a maximum. 



to 




part v 

Fig. 4. Left Panel: Optical eccentricity e pt (solid curve) and 
participant eccentricity ep ar ^ (dash-dot curve) [10] vs. partici- 
pant nucleon number ?ip ar ^. Right Panel: The same curves as 
Fig. |3] vs. mean participant path length v. The vertical dash- 
dot line denotes the mean value of v for N-N collisions |12] 



5 Eccentricity 

Interpretations of the Au-Au azimuth quadrupole require 
understanding the geometry of the initial collision system 
as characterized by the eccentricity |17] defined by 



of this calculation given by 



(yY 



(vY 



(6) 



where x and y are coordinates in the plane perpendicu- 
lar to the beam axis and x refers to the direction along 
the reaction plane. However there are major uncertainties 
about the most appropriate way to calculate eccentricities 
for nuclear collisions. 

In the optical model the eccentricity is approximated 
by a continuous transverse density profile, typically calcu- 
lated from a Woods-Saxon distribution. We shall refer to 
this as e pt- The participant model of eccentricity uses a 
Monte Carlo Glauber calculation that builds nuclei by ran- 
domly placing nucleons according to a Woods-Saxon dis- 
tribution [18]. We shall refer to this as £p ar t- This model is 
favored by many both because of the apparent scaling the 
quantity V2 / £part exhibits between Cu-Cu and Au-Au col- 
lision systems |19j and because it provides a mechanism 
for flow fluctuations. However, in 19] V2 was measured 
with the event plane method which does not adequately 
account for nonflow effects which are large in central col- 
lisions and as stated in section [4] there is reason to believe 
that flow fluctuations may not play a large role. If nonflow 
contributions have not been completely removed then the 
scaling behavior can change when they are. 

Figure |4] shows the differences between the participant 
and optical eccentricities. The optical calculation was done 
according to the Jacobs and Cooper Wounded nucleon 
method with a Woods-Saxon potential [20], however the 
calculation was done using a nucleon-nucleon cross section 
of 42 mb which is more appropriate for RHIC collision 
energies. We found that there is a simple parametrization 



e opt( n b 11 



f 



5.68 



log 
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bin 



logi 



1136 
"bin 



(7) 



The large value of £p ar t for the most peripheral col- 
lisions is unphysical as it implies that the nucleons are 
point-like objects acting at a distance. The nonzero ep ar ^ 
for the most central collisions seems to give desirable scal- 
ing properties when compared to event-plane V2 [E] ■ How- 
ever, our measurement of 1>2{2_D} goes to zero for central 
collisions, which favors e pi- Given the possibility that el- 
liptic flow is a "long-wavelength" probe that is dependent 
only on initial-state collision parameters [10j . we will use 
e opt- 



6 Quadrupole Scaling 

Returning to the quadrupole data in Fig. [2] we note two in- 
teresting features: 1) all energies are described by the same 
centrality variation (dashed curves), and 2) the energy 
dependence of the quadrupole amplitude is propotional 
to \og(y/s n jv/13 GeV). A similar energy dependence was 
observed for (p t ) fluctuations/correlations attributed to 
minijets [7]. 

In Fig.01/e^p t -Z\ j o[2]/ v /p^| is plotted vs. the number 

of binary collisions, nu^ n , times a parameter describing the 
energy scaling 



iZ(V«Aw) = log{Vi^/13 GeV}/ log(200/13). (8) 
For all Au-Au collisions the data can be described by 

^=(%iii> y/sNN) = AR{y/s NN )n hin e 2 t (n bin ) (9) 
V^ref 

where coefficient A is defined by lOOOA = 4.5 ±0.2. Devi- 
ations of 17 GeV event-plane V2 from the linear trend are 
consistent with expected contributions from minijets (10] . 
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Fig. 5. The azimuth quadrupole component divided by optical 
e 2 vs. energy-dependent factor R(y/SNN) times 200 GeV Au- 
Au binary-collision number n^ n (b). The data are consistent 
with simple proportionality over three decades. The dash-dot 
curve is obtained by using £p ar t- 



Equation © accurately describes measured p t -integrated 
azimuth quadrupole moments in heavy ion collisions for 
all centralities down to N-N collisions and all energies 
down to ytspjN ~ 13 GeV. Transformed to each plotting 
space it defines the dashed curves in Fig. [2] The dash-dot 
curves passing through the 200 GeV ID projection points 
in Fig. [3 are obtained by adding "nonflow" parametriza- 
tion 92 /2tt = 0.004I/ 1 - 6 to Ap[2]/ y /p^ [10]. 

We can contrast this simple description of the data 
with the hydro-inspired v^/c vs. 1/ Sdn c h/dn format [2T] 
shown in Fig. [6] For thermal equilibrium the expectation is 
the ideal- hydro limit v-ijt — ► constant. Previous vi mea- 
surements have been interpreted to suggest that central 
Au-Au collisions at 200 GeV attain the ideal hydro limit 
(full thermalization over some substantial volume) [22] . 
The present analysis is inconsistent with such expectations. 
We find a universal linear trend on nuj n and ^/snn for 
all data as in Fig. [5] 



7 Minijets 

We have also fit the same-side minijet peak in the same 62 
and 200 GeV Au-Au autocorrelations. The simple behav- 
ior of the quadrupole component can be contrasted with 
the minijet peak, which exhibits a remarkable transition 
in excess of binary collision scaling around v = 2.5 at 200 
GeV and v = 3.4 at 62 GeV [23] as seen in Fig. [7] and 
Fig. [5] This deviation is a strong indication of medium ef- 
fects. In more peripheral cases the behavior approximates 



co 0.3 
0.25 
0.2 
0.15 1 
0.1 
0.05 
ft 



~i 1 1 1 1 1 1 1 r 



hydro 



200 GeV * 



-♦62GeU 
17 Ge\fe 



STAR 
Preliminary 



LDL 



J I I I I I I I L 







10 



20 



1/S dn h /dr| 



Fig. 6. A hydro-inspired plot of vije vs. LDL (low-density 
limit) parameter 1/S dn c h/dn [21] . 



a simple superposition of nucleon-nucleon collisions. In the 
quadrupole component no such transition is observed. We 
have already seen in Fig. [2] that there is a smooth evo- 
lution according to the collision geometry from the most 
peripheral to the most central collisions, and this evolu- 
tion is described by the parametrization in Eq. @. 

From our fits we conclude that the m = 2 Fourier 
component of the Gaussian minijet peak is the primary 
source of nonflow effects in V2{\D} and similar measures. 
The trend in Fig. [7] also implies that the magnitude of the 
nonflow contribution is strongly centrality-dependent. 



8 Conclusion 

In conclusion, azimuth quadrupole moments vs. centrality 
for Au-Au collisions were measured at 62 and 200 GeV 
separately from minijet contributions (nonflow) via fits 
to 2D angular autocorrelations. These fits use the shape 
of the correlation structures to quantify pair-wise relative 
pseudorapidity dependent nonflow contributions to two- 
particle V2- 

Using the optical Glauber model of eccentricity we find 
that the trend Ap[2]/^p^ oc R(y/s NN )nbi n (b)e^ t (b) 

with i?(- v /sjvjv) oc \og(y/sN at/13 GeV) describes the data 
for all energies and centralities. All p t -integrated Au-Au 
azimuth quadrupole data from 13 to 200 GeV are rep- 
resented by two initial-state parameters. This description 
does not rely on the standard hydrodynamical expectation 
of an equation-of-state, which we see no evidence for. 
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Fig. 7. Minijet peak amplitudes vs. mean participant path 
length v for 200 GeV Au-Au collisions (solid dots) and 62 GeV 
Au-Au collisions (open dots). The dashed and dotted lines rep- 
resent binary collision scaling estimates. 
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Fig. 8. Minijet peak widths in r\A vs. v. 
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